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This article considers dynamical entanglement in non-relativistic particle scattering. Three ques- 
tions are explored: what kinds of entanglement occur in this system, how do global symmetries 
constrain entanglement, and how do the boundary conditions of scattering affect dynamical entan- 
glement? First, a simple model of scattering spin systems is considered, then the full system is 
discussed. 
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I. INTRODUCTION 



In scattering, the initial state of each beam (or the beam and the target) is represented by a density operator 
constructed from eigenvectors of the free Hamiltonian. The total in-state is the direct product of these density 
matrices and is separable. The interaction provides the mechanism for entanglement and the final state is generally 
not separable. This paper studies this kind of dynamical entanglement in scattering. This subject has recently 
received increased interest 0, 0| . 

In particular, this work investigates how entanglement is dynamically generated in elastic, non-relativistic, two- 
T-H particle scattering processes. The following questions will receive attention: (1) What kinds of entanglement occur 
in non-relativistic particle scattering? Non-relativistic particles have continuous (e.g., momentum) and discrete (spin 
component) degrees of freedom, and therefore the state space has a richer structure and allows for different mani- 
festations of entanglement. (2) How do conservation laws limit the dynamics of entanglement? Since the dynamical 
evolution of the system must respect global Galilean invariance, not all entangling operations can arise through scat- 
tering processes. (3) How do the boundary conditions of scattering affect the dynamics of entanglement? The in-state 
must be separable. What other requirements on the interaction and on the in-state are necessary or sufficient so that 
entanglement is generated? 

These questions do not have complete answers. In what follows, techniques from symmetry representations in 
quantum mechanics, like Clebsch-Gordan methods, will be combined with the S-matrix formulation of scattering 
' dynamics to give partial answers and uncover new avenues for further inquiry. Only pure states will be considered. 



II. SCATTERING SPIN SYSTEMS 



^ , To develop intuition and some terminology, we first consider a simpler model of scattering. We will neglect the 
momentum degrees of freedom of non-relativistic particles and consider only the effect of scattering on the spin 
degrees of freedom. In other words, we can consider two spin systems interacting for some finite time interval by a 
rotationally-invariant interaction (the relevant global symmetry). As will be discussed below, low-energy scattering 
by a central interaction approximates the model of scattering spin systems. 

For simplicity, we consider the scattering interaction of two spin-1/2 systems, particle A and particle B. The Hilbert 
space for two spin-1/2 systems is the tensor product of each particle's Hilbert space Since the free 

systems before and after the interaction have rotational symmetry, each is the representation space for the unitary 
irreducible representation (UIR) D^^/^^ of SU(2). For the in-state and out-state, we will use the spin component basis, 
i.e. 

Si\±)A^±ih/2)\±)A. (1) 

We assume that the two systems are initially in pure states of spin and that they have never interacted before. The 
coordinates can always be chosen so that the initial state of particle A is — \+)a- Then, except for an overall 
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phase, the most general state for the second particle is 

WS) = cos(0/2)|+)b + e^^'^ sin(0/2)|-)B, (2) 

where Q < <tt and Q < (j) < 47r. The total in-state is the direct product 

|<^™) = cos(0/2)| + +) + 6^-^/2 sin(0/2)| + -), (3) 

where, for example, | + +) = \+) A ® This is a stationary state of the free Hamiltonian before the scattering. 

Scattering brings the spin degrees of freedom into contact with each other. After the interaction is over, the new out- 
state 10°"*) will again be a stationary state of the free Hamiltonian. The dynamics of scattering can be represented by 
the unitary S-operator which transforms the in-state into the out-state S| = | As a consequence of rotational 

symmetry, the S-operator must satisfy the relation 

[S,S]=0 (4) 

where S is the total spin operator with component operators Si — Sf (g)! + I (g) Sf . Since S is a vector operator 
and the S-matrix is a scalar operator, the Wigner-Eckart theorem tells us that the S-matrix can only be a function of 
total spin s. Therefore, dynamics will be more transparent if we decompose the total Hilbert space into a direct sum 
of UIRs of the global symmetry g) = © 7i*^^. The in-state (EJ is transformed into the total spin basis 

jssa) using the standard Clebsch-Gordan coefficients for the rotation group to get 

= cos(0/2)|ll) + i=e**/2sin(^/2) (|10) + |00)) . (5) 

Unitarity constrains the S-operator and the S-matrix in the direct sum basis has the form 

{s.ss\S\s's',)=e^^'^Sss'S,.,,,, (6) 

where are the phase shifts. All possible dynamics in spin- 1/2 scattering systems can be specified by two parameters 
5o e [0, tt) and 5i e [0, tt). Further, we expect that a global phase is irrelevant and so only the difference {Sq — (5i) G 
(— 7r,7r) should affect physical observations. 

Using and 0, the most general out-state (within the chosen coordinate system) is 

10°"*) = S-l^™) 

= cos(6'/2)e2'*i|ll) + ^e"^/2sin(6'/2) (e^'^^ |10) e^'^^lOO)) 
v2 

= cos(0/2)e2'^n + +> (7) 

-hie*'*/^ sin(6'/2) {(e^'^i + e^'^°) \ + -) + {e^'^' ~ e^'^«) \ ~ +)} . 

So for this simple model, four parameters {0, 9, i5o, i5i} are sufficient to describe every possible in-state (the first two 
parameters) and all dynamics arising from central forces (the second two). We now turn to the three questions 
mentioned in the introduction. 

What kinds of entanglement occur in non-rclativistic particle scattering? In this model, the global Hilbert space 
is finite dimensional and can be factored into the direct product of two Hilbert subspaces. Therefore, the entropy of 
entanglement (EoE) provides an unambiguous measure of entanglement in pure states. For the full case of momentum 
and spin degrees of freedom, the answer to this question is more complicated, as will be discussed below. 

Calculating the EoE of Q, we find 

i?(|0™*)) - 1 - log2 [(1 + x)(i+-)(l - x)(i--)] , (8) 

where 



X = \/l - sin'^(6'/2) sin2(2A(5) (9) 

and AS = 5q — di. The EoE (jSJ has a maximum E ~ 1 when a; = and a minimum E = (no entanglement) when 
X = 1. 

How do symmetries limit the dynamics of entanglement? In this case, rotational symmetry means that the S-matrix 
must have the highly- restricted form ©. From ||SJ), we see that the phase shift difference AS determines whether 
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the dynamics can be entangling. In particular, we see that only for Ai5 = ±7r/4 can the S-operator lead to maximal 
entanglement. Such operators have been called perfect entanglers |14| because they are operators which can take an 
unentangled vector to a maximally-entangled vector. 

How do the boundary conditions of scattering affect the dynamics of entanglement? The EoE of the in-state is zero 
by construction as a consequence of the assumption of inital independence of the particles. For which initial states are 
there perfectly-entangling dynamics? In other words, for which values of (f) and 9 are there solutions to E{\(t)°'^*)) = 1? 
The phase (j) is unimportant, and only when = tt is the in-state = | H — ) maximally-entangleable. This in-state 
has equal components in both partial waves. Even then, maximal entanglement can only occur for A(5 = ±7r/4. In 
this case, the out-state is 

|,/>°"*)=c.i=(|+-)T *!-+)) (10) 

where there is an overall phase lo = exp(i0/2) exp(2i(Si) exp(±i7r/4). 

The lessons of the toy model can be summarized as follows: Enforcing global symmetry limits the set of perfectly- 
entangling dynamical operators and the sets of separable in-states that can be perfectly entangled. For some in- 
states, no entanglement is possible for any dynamics. Finally, only a limited set of maximally- entangled out-states 
are accessible via dynamical entanglement. 



III. THE FULL PROBLEM 



Next we discuss the problem of entanglement in non-relativistic scattering. We will attempt to show that all of 
the above observations also appear to be true for the full problem, although there are many complications and many 
other kinds of entanglement phenomena can occur. 

The first step is to cast the mathematics in terms of the language of symmetry. The symmetry group is now the 
Galilean group Q (or, more properly, its quantum-mechanical covering group Q). UIRs of Q are associated with a 
single, free, non-relativistic particles They are denoted and are characterized by invariant mass to, 

internal energy w, and spin s. The representation space is denoted T^^'"''"'*). Eigenkets of momentum and spin 
component Ip, S3 (toius)) G (<i)(™'"'''*))^ form a basis for a dense subspace <I)(™''"'*) c 7i('"''"'*) of well-behaved vectors 
(see Ref. (ual) for details). Unlike in the case of relativistic particles, the space $(™''"'*) can be factorized into the 
direct product of momentum and spin degrees of freedom: 

<j,(m,to,s) ^ (^(m,w) ^ j^{2s+l)^ ^^^^^ 

It is assumed that a pure in-state is separable, i.e. |(/)™) = |</>™) (X" 10^) where |(/)™) £ <j)(™i^«'i^si)_ rpj^g S-operator 
transforms to |(/)°"*), and if the scattering is elastic, then we have 

l^out^ g <j,(mi,t"l,si) ^ <j,(m2,t"2,S2)^ (^12) 

but in general it will no longer be separable. Because of global Galilean invariance, the S-matrix will have a very 
restricted set of non-zero matrix elements, just as in the spin system model. To exploit this, we can use Clebsch- 
Gordan methods to separate the direct product of single-particle UIRs into a direct sum of UIRs labeled by the global 
invariants: 

^{mi,wi,si) ^ ^{m2,W2,S2) _ <j,(M,Wj) [r;i ,772] ^ ^-j^g^ 

where the set {M, W,j} are the labels of the UIR and are global Galilean invariants with ranges that depend on the 
sinle particle invariants 15]. The same set of invariants may appear more than once in the direct sum; rji, r]2 label the 
degeneracy. In standard analysis, these are chosen to be the total spin angular momentum s and the orbital angular 
momentum I and this process is called partial wave analysis. 

The Clebsch-Gordan coefficients can be used to transform the wave functions between the direct product basis 
, s'^{mAWASA)) , sf {rriBWBSB)) and the UIR direct sum basis with spin-orbit coupling \pj3{MWj)[ls]) = 

\pj3n). In the direct sum basis, because of Galilean invariance, the S-matrix will have the form .17\ . 

(pj3«|S|p'j>') = NiM,W)6'ip - p')6{W - W')S,,,S,,,,{ls\\S^'^^^ Ill's'), (14) 

where N{M,W) is determined by the normalization of the basis kets and Sj^j'^{ls\\S''^^^\\l's') is the reduced matrix 
element that mixes spin and orbital angular momentum. 
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What kinds of entanglement occur in this kind of scattering? A much wider variety, because when two non- 
relativistic particles interact, the total Hilbert space cab be decomposed into the direct product of four independent 
subspaces, two for each particle's spin and two for momentum. Tracing out the degrees of freedom of only one of 
these constituant spaces could result in a mixed state, pointing toward entanglement. Taking a trace over different 
combinations could lead to EoEs associated with two-particle spin entanglement, two-particle momentum entangle- 
ment, internal entanglement between one particle's spin and momentum, and mixed varieties. Even in-states with 
independent particles could have entanglement present in the form of internal spin- momentum entanglement. 

Just like the spin system model, where global rotational symmetry limits when and how dynamical entanglement 
occurs, global Galilean symmetry has consequences for the dynamics and entanglement: not all unitary transforma- 
tions are possible, so not all interactions are entangling. Further, not all separable in-states may be entangleable. 
The interplay between the symmetry of the dynamics and the boundary conditions of scattering could lead to some 
states that do not become entangled by any interaction, others that entangle under a wide class of interactions, and 
exhanges between different kinds of entanglement. 

We conclude with one simple example that connects the full problem back to the spin systems model. In the case of 
central interactions, I and s are also separately conserved. Then a more convenient basis to understand the dynamics 
of entanglement is the l-s basis |pZ3S3(MT/F)[^s]) = \pl3S3n), and the S-matrix elements in this basis are ^3 

(pZ3S3n|S|p';^4n') = N{M, W)6'ip - p')6{W - W')6u'Si,i>Sss'S,,,'^e^''^'^^'' . (15) 

If the in-state is sharply peaked around a particular total energy W, and only one orbital angular momentum par- 
tial wave contributes to the scattering appreciably (as in the low energy limit), then the spin-dependent part of 
the interaction can be separated from the momentum-dependent part and the toy model above becomes a realistic 
approximation 

However, this is only the simplest case. Many new phenomena and types of entanglement wait to be explored. 
There are complications, namely working with continuous degrees of freedom, but recent progress with Gaussian 
states makes these difficulties appear tractable. 
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